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We investigate the ground state phase diagram of the half-filled t — t' repulsive Hubbard model in 
the presence of a staggered ionic potential A, using the continuum-limit bosonization approach. We 
find, that with increasing on-site-repulsion U, depending on the value of the next-nearest-hopping 
amplitude t' , the model shows three different versions of the ground state phase diagram. For 
t' < t't, the ground state phase diagram consists of the following three insulating phases; Band- 
Insulator at [/ < Uc, Ferroelectric Insulator ai Uc < U < Us and correlated (Mott) Insulator at 
U > Uc- For t' > t'c there is only one transition from a spin gapped metallic phase ai U < Uc 
to a ferroelectric insulator at f/ > f/c. Finally, for intermediate values of the next-nearest-hopping 
amplitude ti < t' < t'^. we find that with increasing on-site repulsion, at Ud the model undergoes a 
second-order commensurate-incommensurate type transition from a band insulator into a metallic 
state and at larger Uc2 there is a Kosterlitz-Thouless type transition from a metal into a ferroelectric 
insulator. 

PACS numbers: 71.10.Fd, 71.27.+a, 71.30.+h 



I. INTRODUCTION 

During the last decades, the Mott metal-insulator tran- 
sition has been a subject of great interest jii^i^ In the 
canonical model for this transition - the single-band Hub- 
bard model - the origin of the insulating behavior is the 
on-site Coulomb repulsion between electrons. For an av- 
erage density of one electron per site, the transition from 
the metallic to the insulating phase is expected to oc- 
cur with increasing on-site repulsion when the electron- 
electron interaction strength U exceeds a critical value 
Uc, which is usually of the order of the delocalization 
energy. Although the underlying mechanism driving the 
Mott transition is by now well understood, many ques- 
tions remain open, especially about the region close to 
the transition point where perturbative approaches fail 
to provide reliable answers. 

The situation is more fortunate in one dimension, 
where non-pcrturbative analytical methods together with 
well-controlled numerical approaches allow to obtain an 
almost complete description of the Hubbard model and 
its dynamical properties*^ However, even in one dimen- 
sion, apart from the exactly solvable cases, a full treat- 
ment of the fundamental issues related to the Mott tran- 
sition still constitutes a hard and challenging problem. 

Intensive recent activity is focused on studies of the 
extended versions of the Hubbard model which display, 
with increasing Coulomb repulsion, a transition from 
a band-insulator (BI) into the correlated (Mott) insu- 
lator phase. Various models considered include those, 
which show a continuous evolution from a BI into the MI 
phase^ii as well as those, where the transformation of 
a BI to a correlated (Mott) insulator takes place via a 
sequence of quantum phase transitionSi ^'^'-'^^d-'^'-'^^d^'-'^'^d^ 



Intensive recent activity is focused on studies of the ex- 
tended versions of the Hubbard model with alternating 
on-site energies ±A, known as the ionic Hubbard model 
(IHM)4^a0iiia2i3^ii5a6i7^i9^ The model has a 

long-term history;^ however the increased current inter- 
est widely comes from the possibility to describe the in- 
teraction driven BI to MI transition within one model. 
In one dimension this evolution with increasing on-site 
repulsion is characterized by two quantum phase transi- 
tions: first a (charge) transition from the BI to a ferro- 
electric insulator (FI) and second, with further increased 
repulsion, a (spin) transition from the FI to a corre- 
lated Ml4 Detailed numerical studies of the ID IHM 
clearly show, that an unconventional metallic phase is 
realized in the ground state of the model only at the 
charge transition point and the BI and MI phases are 
separated in the phase diagram by the insulating fer- 
roelectric phasCi ^d^d^d^ Studies of the 2D IHM using 
the cluster dynamical mean field theory, reveal a similar 
phase diagram 

On the other hand, recent studies of the IHM us- 
ing the dynamical mean field theory (DMFT) approach, 
show that in high dimensions the BI phase can be sepa- 
rated from the MI phase by the finite stripe of a metallic 
phascii^ Moreover, recent studies of the IHM with site 
diagonal disorder using the DMFT approach, also show 
the existence of a metallic phase which separates the BI 
phase from the MI phase in the ground state phase dia- 
gram of the disordered IHMi^i It looks so, that in low- 
dimensional models with perfect nesting of the Fermi sur- 
face, the metallic phase is reduced to the charge transi- 
tion line, while in higher dimensions the space for real- 
ization of a metallic phase opens. Note that the very 
presence of a metallic phase along critical lines separat- 
ing two insulating phases is common for ID systems with 
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competing short-range interactions responsible for the 
dynamical generation of a charge gapi^ However, one- 
dimensional models of correlated electrons, showing with 
increasing electron-electron coupling a transition from an 
insulating to a metallic phase are less knowni^ 

In this paper we show, that the one-dimensional 
half-filled ionic-Hubbard model supplemented with the 
next-ncarcst-ncighbor hopping term {t') is possibly the 
simplest one-dimensional model of correlated-electrons, 
which shows a ferroelectric ground state in a wide area 
of the phase diagram easily controlled by the model pa- 
rameters. 

We also show that in a certain range of the model 
parameters the t — t' ionic-Hubbard chain shows, with 
increasing Coulomb repulsion, a transition from a band- 
insulator to a metal and, with further increase of the 
Hubbard repulsion, a transition from a metal to a ferro- 
electric insulator (see Fig.l). 




Uci U 

FIG. 1: Qualitative phase diagram of the half-filled t — t' 
ionic-Hubbard chain for t' > 0.5t and weak and moderate 
values of t he on-site Hu bbard repulsion U. The parameter 
t'^ = O.bty^l + {A/4ty + A/8 corresponds to the insulator 
to metal transition point in the free ionic chain, where A/2 
is the amplitude of alternating ionic potential. The dashed 
curve line marks the I-M transition line in the t — t' Hubbard 
model, corresponding to the case A = 0. The dashed line 
t' — Q.5t is given as an eye guide. 



The paper is organized as follows. In Section II, the 
model and its several important limiting cases are briefly 
reviewed. In the Section III the weak-coupling bosoniza- 
tion description is obtained. The results arc summarized 
in Section IV. 



II. THE MODEL 

The Hamiltonian we consider is given by 

n,(T 

+ Y.y'^'+ y) '^"■'^ + ^ 5Z /^".TPn,i . (1) 

Here c\^^ (c„_o-) ^^^^ electron creation (annihilation) op- 
erators on site n and, with spin projection a =tj-li 
Pn.cr — clia^^ncr- ^hc ncarcst-ncighbor hopping ampli- 
tude is denoted by t, the next-nearest-neighbor hop- 
ping amplitude by t' {t,t' > 0), A is the potential en- 
ergy difference between neighboring sites, U is the on- 
site Coulomb repulsion and the band-filling is controlled 
by the proper shift of the chemical potential 5^. For 
i' = 0, we recover the Hamiltonian of the ordinary ionic- 
Hubbard chain, while for A = the Hamiltonian of the 
t — t' Hubbard chain. 

At [/ = the model is easily diagonalized in momen- 
tum space to give a dispersion relation in the first Bril- 
louin zone — 7r/2 < fc < 7r/2 

E±{k) = 2t'cos2fc- (5Ai± v/4t2 cos^ fc + (A/2)2 . (2) 

Let us first analyze the dispersion relation ([2]). For 
t,t' > the absolute maximum of the upper band is 
reached at fc = 0, while the absolute minimum 

i?+„ = -2t'-5M + A/2, (3) 

at fc = ±7r/2. 

The lower band shows a more complicated dependence 
on the model parameters. For 

t' < t'^ = OM v/l + (A/4i)2 - A/8 (4) 

the absolute minimum of the lower band is reached at 
fc = and while the absolute maximum at fc = ±7r/2 
and is equal to 

E- ^ ^2t' - 5fi - A/2 . (5) 

In the case of half-filling the lower band is completely 
filled and the upper band is empty. The system is a 
band-insulator with a gap in the excitation spectrum 

Ae.c = - i?„-a. = A. (6) 

The corresponding shift of the chemical potential Sfi = 
—2t' is easily determined from the condition i?^j„ -I- 
^max — 0- Thus, for t' < t'^ the ground state and low- 
energy excitation spectrum of the model are not affected 
by the increase of the next-nearest hopping amplitude t' . 
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The effect of t' becomes nontrivial for t' > t'^, when 
the absolute maximum of the lower band 



E-^,, = 2t' - <5m - v/4t2 + (A/2)2 



is reached at fc = 0. For 



t' < t'^= 0.5i ^Jl + {/^/AtY + A/8 



(7) 



the absolute maximum of the lower band at fc = re- 
mains lower than the absolute minimum of the upper 
band at fc = 7r/2. Therefore the system remains an insu- 
lator, but the indirect gap 



A 



As = - + 



(A/2)2 -4t' 



(8) 



decays linearly with increasing t' and finally vanishes at 
t' = t'^. It is straightforward to find, that the corre- 
sponding shift of the chemical potential which ensures 
half-filling in this case is given by Sfi = —2t'^. 

At t' = t'^ the t~t' ionic chain experiences a transition 
from a band-insulator to a metal. For t' > t'^ a gap- 
less phase is realized, corresponding to a metallic state 
with four Fermi points ±fc^i and ±kp2, which satisfy the 
relation 2(fc;'i — fci?2) = tt. 

At J7 > several limiting cases of the model ([ij have 
been the subject of intensive current studies. In particu- 
lar intensive recent activity has been focused on studies 
of the ground state phase diagram of the ionic-Hubbard 
model (IHM), corresponding to the limiting case t' — 
and of the t — t' Hubbard model corresponding to the 
limiting case A = 0. 

Current interest in this model mostly originates 
from the possibility to describe the evolution from the 
band insulator (BI) at [/ ^ A into a correlated 
(Mott) insulator (MI) ioi U ^ A within a single 
system^iSiiadi^aSiiiasae^as^ In the case of the one- 
dimensional ionic-Hubbard model this evolution is char- 
acterized by two quantum phase transitions in the ground 
state<^ With increasing U the first is a charge transition 
at J7 = Uc, from the BI to a bond-ordered, spontaneously 
dimerized, ferroelectric insulator (FI). At the transition 
point the charge sector of the model becomes gapless, 
however for U > Uc the charge gap opens again. When 
U is further increased, the second (spin) transition from 
the FI phase into the MI phase takes place at Us > Uc- At 
this transition the spin gap vanishes and the spin sector 
remains gapless in the MI phase for U > Ug. A similar 
ground state phase diagram has been recently established 
also in the case of 2D IHM using the cluster dynamical 
mean field theoryji^ Thus, for low-dimensional versions 
of the IHM with perfect nesting property of the Fermi 
surface, the transition from a BI to a FI is characterized 
by the presence of a metallic (charge gapless) phase only 
at the transition line. In the ground state phase diagram 
of the BI phase is separated from the MI phase by a FI 
phase.^'ia^^ 

The homogeneous half-filled t — t' Hubbard 
chain is a prototype model to study the metal- 
insulator transition in one dimension and therefore 



has been the subject of intensive studies in recent 
yearsffi^iM^i^^^^ia^^ As in the case of the 
ionic-chain, in the t — t' Hubbard model an increase 
of the next-nearest hopping t' changes the topology of 
the Fermi surface: at half-filling and for t' < 0.5i, the 
electron band of has two Fermi points at kp ~ ±tt/2, 
separated from each other by the umklapp vector g = tt. 
In this case, a weak-coupling renormalization group 
analysis predicts the same behavior as for t' = - the 
dynamical generation of a charge gap for U > 0, and 
gapless magnetic excitations <^ 

For t' > Q.5t, the Fermi level intersects the one-electron 
band at four points (fc^ ^ ±7r/2). For weak Hubbard 
coupling ([/ <^ t) the infrared behavior is governed by the 
low-energy excitations in the vicinity of the four Fermi 
points, in full analogy with the two-leg Hubbard model4^ 
The Fermi vectors kp are suflticiently far from tt/2 to 
suppress first-order umklapp processes and therefore the 
charge excitations are gapless, however the spin degrees 
of freedom are becoming gapped i^'^^'^i^S'^ Since at 
half-filling 4(fc^ ~ kp) = 27r, with increasing on-site 
repulsion higher-order umklapp processes become rele- 
vant for intermediate values of U. Therefore, starting 
from a metallic region for small [/ at a given value of t' 
{t' > 0.5t), one reaches a transition line U = Uc{t') above 
which the system is insulating with both charge and spin 
gapsi^ 

As we will show in this paper for t' < t'^, where 
the topology of Fermi surface is restricted to two Fermi 
points, the ground state phase diagram of the t — t' ionic- 
Hubbard chain coincides with that of the IHM, while 
for t' > <J,, where the model is characterized by four 
Fermi points - it coincides with that of the t — t' Hub- 
bard chain. Most interesting is the case t'^ < t' < t'^, 
where with increasing on-site repulsion we observe two 
transitions in the ground state: at Ud the model under- 
goes a second-order commensurate- incommensurate type 
transition from a band insulator into a metallic state and 
at large Uc2 there is a Kosterlitz-Thouless type transition 
from a metal into a correlated ferroelectric insulator. 



III. BOSONIZATION RESULTS 

In this section we analyze the low-energy properties 
of the t ~ t' ionic-Hubbard chain using the continuum- 
limit bosonization approach. Wc first consider the regime 
U, A, t' <^ t, linearize the spectrum in the vicinity of the 
two Fermi points kp = ±tt/2 and go to the continuum 
limit by substituting 



(9) 



where x = tioq, ao is the lattice spacing, and ipRa-ix) and 
'>pLa{x) describe right-moving and left-moving particles, 
respectively. The chosen type of decoupling of the model 
into " free" and " interaction" parts allows to treat the gap 
"creating" (A and U) and gap " dcstructing" {t') terms on 
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equal footing and reveals easily their competition within 
the continuum-limit treatment. 

The right and left fermionic fields are bosonizcd in the 
standard way4^ 

1 



„JV^[$„(x)+e„(x)] 



1 

I 



(10) 
(11) 



where ^a{x) and Qaix) are dual bosonic fields, dt^a = 
vpdxQcr and vp = 2tao- 

This gives the following bosonized Hamiltonian: 



where 
and 



-^a,$.-7r^sinV4^$4 ('T=T,i)(12) 

H — ^ sin V47r<l>f sin V47r$x ] |. (13) 
Here we have introduced 

for t' < 

2 (f - <) for it <t' <t[ ■ 

The next step is to introduce the charge 

^c^^iH+cPl), = ^(01-01) (15) 

and spin fields 

vs = ^{^i-^i), iDs = ^{9^-eo (16) 

to describe corresponding degrees of freedom. After a 
simple rescaling, we arrive at the bosonized version of 
the Hamiltonian ([T|) 

TL — TLs + TLc + , 

where 



(14) 



Us 



dx 



fJ-cSX -dxipc - o ^ 2 coa{V8nifc) \ , (17) 

V TT ZTT^Oq J 

dx{|[(a.^.)2 + i(a,,?,)2] 

-^cos(V8^^,)|, (18) 



Tics = J dx sin (^\/27r(fc^ cos ^"\/27r(^s^ (19) 

with the charge stiffness parameter Kc ~ 1 — C// Ant at 
t//47rt < 1. 



A. Non-interacting case 

To assess the accuracy of the continuum-limit treat- 
ment it is instructive to start our bosonization analysis 
from the exactly solvable case of the ionic-chain. 

At [/ — the system is decoupled into the "up" 
and "down" spin component parts Ti. = + H^, 
where for each spin component the Hamiltonian is the 
sine-Gordon model with topological term Each 
of these Hamiltonians is the standard Hamiltonian for 
the commensurate-incommensurate transition, which has 
been intensively studied in the past using bosonization^ 
and the Bethe ansatz4i This allows to apply the theory of 
commensurate-incommensurate transitions to the metal- 
insulator transition in the considered case of a half-filled 
t — t' chain with ionic distortion. 

At fj.eff = 0, the model is described by the the- 
ory of two commuting sine-Gordon fields (sin/3$o-) with 
= Att. In this case the excitation spectrum is gapped 
and the excitation gap is given by the mass of the " up" 
("down") field sohton M| = = A/2. In the ground 
state the $| and $| fields are pinned with vacuum expec- 
tation values (0|$ct|0) = -\/7r('^+l/4)- Using the standard 
bosonized expression for the 2kp modulated part of the 
charge density^ 

Pc{x) ~ (-l)"^^sin(V4^<i>,(z)) 



TTflo 



(-1)" 

27rao 



sin(v27riy9c(a;)) cos(v27r(ps(a;)) (20) 



we obtain that at fi^f / = the ground state of the system 
corresponds to a CDW type band-insulator with a single 
energy scale given by the ionic potential A. 

At fJ-eff 7^ it is necessary to consider the ground 
state of the sine-Gordon model in sectors with nonzero 
topological charge. The competition between the chemi- 
cal potential term (t') and the commensurability energy 
given by A finally drives a continuous phase transition 
from a gapped (insulating) phase at p^jj < ^'^^g to a 
gapless (metallic) phase at 



(21) 



Using (|T4l) we easily obtain, that the critical value of the 
n-n-n hopping amplitude t' , obtained from the condition 
coincides with the exact value for the ionic-chain 
given in ([7]). 

As we observe, the insulator-metal transition at t' > t'^ 
is connected with a change of the topology of the Fermi 
surface and a corresponding redistribution of the elec- 
trons from the lower ("-") band into the upper ("+") 
band. We use as an order parameter of this transition 
the number N+ of electrons transferred into the "-I-" 
band, which is related to the value of the new Fermi 
point /cq ^ \/ ~ Mc- At the transition point the com- 
pressibility of the system is 



(22) 
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showing an inverse square-root singularity, where Eq is 
the ground state energy. 

Before we start to consider the interacting case, it is 
useful to continue our analysis of the noninteracting case, 
but within the basis of the charge and spin Bose fields, 
which is more convenient for interacting electrons. 

The Hamiltonian we have to consider now is given by 



Here 



A 

nao 



sin (^\/2TTipc^ cos (^V^nifis^ | • (23) 



We decouple the interaction term in a mean-field man- 
ner by introducing 



(24) 
(25) 



nic — A • (cosv27r(/3s 
rus = A • (sin \p2/Kipc 



and get the mean-field bosonized version of the Hamilto- 
nian 7i = Tic + 'yis which is given by the two commuting 
quantum sine-Gordon models 



Us 



rfa;{^P,,?,)2 + (a,^,)2] 
fJ-cff\ -dxfc - sin(y27rv'c) } , (26) 

V TT TTflo J 

cos(-\/27r(^s)| . 

TTflo J 



(27) 



Although the mean-field Hamiltonian is once again 
given by the sum of two decoupled sine-Gordon models 
(see Eq. (HU), the dimensionality of the cos(/3(/)) oper- 
ators at 0^ = 27r and = 47r are different. Therefore, 
in marked contrast with the bosonized theory in terms 
of "up" and "down" fields p2)) . the pair of Hamiltonians 
given by (|26 p -(|27 p represents a very complex basis to de- 
scribe the BI phase i.e. the CDW state with equal charge 
and spin gaps Ac = A^ = A/2. 

Nevertheless, below we use the benefit of the exact so- 
lution of the sine-Gordon model to get a qualitatively and 
almost quantitatively accurate description of the problem 
even in the "spin-charge" basis. To see this, let us start 
from the case jioS = 0. We will use the following exact re- 
lations between the bare mass m and the soliton physical 
mass M for the sine-Gordon theory with 0^ = 2tt— 



M/A = Co (to/A) 



2/3 



(28) 



and the exact expression for the expectation value of the 
cos 06 field^ 



(cosV2^(p) =Ci (il//A)^/^ . 



Co 



2r(i/6) 
0?r(2/3) 



r(3/4) 
2r(i/4) 



(30) 



and 



Ci = 



2 fSirV^'^ r(3/4) 

rW4) 



(31) 



and A = 2t is the bandwidth. 

Using (I28l)-(13T]) one easily finds that 

Ac/A ^ CoiA/kf'^ {cosV2^ip,)'^/^ 
= CoCf/^(A/A)'/^(M,/A)'/^ 

A,/A = Co(A/A)2/^(sinV2^^c)^/^ 
= CoCi/^(A/A)'/^(A/c/A)'/^ . 



(32) 



(33) 



The self-consistent solution of the Eqs. ([32 |) -([33 |) gives 



Ac = A, = 7A/2 



(34) 



0.94256. 



Thus the mean-field 



(29) 



with 7 = 2Cy^Ci 
treatment of the ionic-chain within the spin-charge basis 
gives not only a qualitatively correct but, rather quan- 
titatively accurate description for the system in the BI 
phase. 

For completeness of our description, let us now con- 
sider the insulator to metal transition in the ionic-chain 
in the charge-spin Bose field basis. The correspond- 
ing mean-field decoupled charge Hamiltonian is given 
by the 0^ = 27r quantum sine-Gordon model with the 
topological term ([26|l . From the exact solution of the 
SG model^ it is known that the excitation spectrum of 
the model at 0^ = 27r consists of solitons and antisoli- 
tons with mass Mc, and soliton-antisoliton bound states 
("breathers") with masses Mc"^^ = 2A/c sin(7r/6) = A/c 
and Afc"=^ = 2Afc sin(7r/3) = VSA/c- The transition 
from the BI into the metallic phase takes place when 
the effective chemical potential exceeds the mass of the 
lowest breather i.e. at /Xe// = Mc = A/2. For /Xe// > Mc 
the vacuum average of the charge field is not pinned at 
all, such that the spin gap {Mg) ~ (sin(-\/27r</'c)) = 0. 

Thus, using the bosonization treatment we easily and 
with good numerical accuracy describe the exactly solv- 
able case of the band-insulator to metal transition, which 
takes place in the t — t' ionic chain with increasing next- 
nearest- hopping amplitude t' . 



B. Interacting case 

At (7 ^ the use a similar mean-field decoupling allows 
us to rewrite the Hamiltonian as two commuting double 
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sinc-Gordon models Ti. = Tic + Ti-s where 

V TT Trao 
-:r^cos(2y2^^,)}, (35) 

YY^^ 1^—^— ]\'_f ^ 

^ cos(V27r(ps) + % cos( VStti^s) f ■ (36) 

Here 

mj; = Ar • (cos-\/27r(ps) , (37) 

= Ar ■ (sill V2Tripc) , (38) 

Mc and Mg are the effective model parameters. The 
renormalized band gap 

A, = A [!-([//[/*)], (39) 

includes the Hartree renormalization of the band-gap by 
the on-site repulsion [/, where U*{A) is a phenomenolog- 
ical parameter. For given A, U* is of the order of Uc in 
the standard IHM with the same amplitude of the ionic 
distortion. The charge stiffness parameter is Kc < 1 in 
the case of repulsive interaction. 

As we observe at J7 > the charge sector is described 
by the double frequency sine-Gordon model with strongly 
relevant basic and marginally relevant double-field oper- 
ators complemented by the topological term. The spin 
sector is also given by the double frequency sine-Gordon 
theory with strongly relevant basic and, at weak-coupling 
([/ ^ 2t), marginally irrelevant double-field operators. 

At ^off = 0, where the charge excitation spectrum 
is gapped, the peculiarity of the charge sector is dis- 
played in the internal competition of the vacuum config- 
urations of the ordered field driven by the two sources 
of gap formation: the ionic term prefers to fix the 
charge field at ^/2ttKc{0\(Pc\0) = tt/2 + 27rn, which cor- 
responds to the maximum of the double-field operators 
~ —U cos{2y/2TrK2ipc) , i.e., it corresponds to a configu- 
ration which is strongly unfavored by the onset of cor- 
relations. On the other hand, the vacuum expectation 
value of the field (0|(/3c|0) = 2?™, which minimizes the 
contribution of the double-field operator for {/ > 0, leads 
to the complete destruction of the CDW pattern, which 
was favored by the alternating ionic potential. 

This type of competition in the double frequency sine- 
Gordon model results in a quantum phase transition from 
the regime where the field is pinned in the vacuum of 
the basic field potential into the regime where the field 
is pinned in the vacuum of the double-frequency cosine 
termi^ Qualitatively the transition point can be esti- 
mated from dimensional arguments based on equating 
physical masses produced by the two cosine terms. This 



allows to distinguish two qualitatively different sectors of 
the phase diagram corresponding respectively to the case 
of weak repulsion (U « A,t), where the ground state 
properties of the system are determined by the band gap, 
and to the case of strong repulsion {U » A, i), where 
the charge sector is characterized by the Mott-Hubbard 
gap. However, the detailed analysis of the critical area in 
the case of the IHM shows, that the BI is separated from 
the Mott insulator by a ferroelectric insulating phase. 




FIG. 2: Qualitative sketch of behavior of the single parti- 
cle (solid line), spin (dashed line) and optical (dashed-dotted 
line) gap as a function of the on-site repulsion U based on the 
exact numerical results obtained in Ref. and Ref. [l3 

One important tool to characterize the different phases 
of the IHM is to study gaps to excited states, in particular 
making contact with the gaps obtained in the bosoniza- 
tion description. Following Ref. [l^ we define the optical 
gap as the gap to the first excited state in the sector with 
the same particle number TV and with Sz = 0, where Sz 
is the z-component of the total spin. The single particle 
gap is determined as the difference in chemical poten- 
tial for adding and subtracting one particle. Finally the 
spin gap is defined as the energy difference between the 
ground state and the lowest lying energy eigenstate in 
the S — 1 subspace. 

In Fig. 2 we present a qualitative sketch of the behav- 
ior of the single particle charge (solid line), spin (dashed 
line) and optical (dashed-dotted line) gap as a function 
of the on-site repulsion U based on the exact numerical 
results obtained in Ref. and Ref. [T^. Two different sec- 
tors of the phase diagram corresponding to the BI and 
MI are clearly shown. These sectors are distinguished 
by the pronounced difference in the U dependence of the 
charge excitation (single-particle excitation) gap. 

Below we treat the ground state phase diagrams of 
the IHM and of the noninteracting ionic-chain as border 
lines of our model. We explore the different character of 
the excitation gap renormalization by the Hubbard re- 
pulsion U and consider the ground state phase diagram 
of the half-filled repulsive t — t' ionic-Hubbard chain. The 
key component of our analysis is based on the assump- 
tion that for arbitrary U the infrared properties of the 
model arc determined by the relation between the con- 
trolled effective chemical potential /ie// and the value 
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of the charge excitation gap. Moreover, in each charge 
gapped sector of the phase diagram, there is only one en- 
ergy scale which is given either by the ionic distortion or 
by the Hubbard repulsion. This allows us to get the qual- 
itative ground state phase diagram, which is summarized 
in Fig. 3. 

At t' < t'^ (i.e. //off = 0) the ground state phase dia- 
gram remains qualitatively the same as at t' = 0: the BI 
phase is separated from a MI phase via a narrow stripe 
of the FI phase. 

At /ieff 7^ 0; but t' < the BI phase is present for 
weak repulsion. With increasing J7 the band-gap reduces, 
which simultaneously leads to a renormalization of the 
effective chemical potential //'' = 2(t' — ij,), where t[. is 
given by ^ with A = A,,. Here we have to consider two 
cases separately. 

At t'^ <t'< 0.5t, the renormalizcd value of the chem- 
ical potential always remains less then the renormalizcd 
single-particle gap A,./2. Therefore with increasing U 
the phase diagram is qualitatively the same as at t' < t'^ , 
i.e. with increasing U the system undergoes a transition 
into a FI phase and with further increase of U into a MI 
phase. Since the effective single band gap for t'^ < t' is 
smaller than the single-particle gap at = 0, the transi- 
tion from a BI into the FI insulator takes place for smaller 
[/, which manifests itself in an extension of the FI phase 
in the direction of lower U . 

At 0.5t < t' < t'^ the BI phase is still realized at 
weak J7, but with increasing repulsive interaction one 
reaches the critical point ?7ci, where the renormalizcd 
value of the chemical potential exceeds the renormal- 
izcd single-particle gap Ar/2. Neglecting quadratic cor- 
rections in A/i, using ((39|) we easily find that Uci — 
U*[l - %{t' - 0.5t)/A] < f7* . At [/ > Uci the BI phase 
is destroyed, the amplitude of the basic frequency field 
operator sm[^/2T:{pc) in the charge Hamiltonian (j35p van- 
ishes and the charge sector of the model is given by the 
Hamiltonian 



Tic ~ I dx 



^— - cos(2V27ri^,^,)| . (40) 



The Hamiltonian (|40p is a Hamiltonian which describes 
the charge sector of the t — t' Hubbard chain at t' > 0.5t 
and contains two regimes of behavior)^ 

a) if the effective chemical potential is larger than 
the value of the correlated (Mott) gap at the transi- 
tion point /icff > A'lciUci) then the BI phase undergoes 
a transition into a charge gapless metallic phase. Since 
in the BI phase the energy scale of the model is given 
by the renormalizcd band-gap, in analogy with the non- 
interacting case of ionic chain we expect that the tran- 
sition from a BI to a metal belongs to the universality 
class of commensurate-incommensurate transitions. 

b) if fj.cS < Mc{Uci) then the BI phase undergoes a 
transition into a charge gapped ferroelectric phase. 



We estimate the charge gap for U ^ t as Mc « 
^r^p-2^t/u g^jjj as Mc w [/ for [/ > i. One finds that for 
t' < t'^ the effective chemical potential is larger than the 
exponentially small Hubbard gap Mc{Uci) and therefore 
for 0.5i < t' < t'^ there opens a window for a transition 
from the BI to a metallic phase with increasing on-site 
repulsion. With further increase of the on-site repulsion, 
at Uc2, when Mc(Uc2) = A*off a- charge gap opens once 
again, and for U > Uc2 the system is in the insulating 
ferroelectric phase (see Fig. 3). 



Ferroelectric 
insulator 




Uc Us 



FIG. 3: Qualitative sketch of the ground state phase diagram 
of tlie t — t' ionic-Hubbard chain in the case of repulsive in- 
teraction. Solid lines mark the phase transitions. 

At t' > t'^ the phase diagram is more simple. Already 
for U ~ Q the ground state corresponds to a metallic 
state, since the effective chemical potential is larger than 
the band-gap. With increasing Hubbard repulsion a tran- 
sition into an insulating phase takes place, when the cor- 
related gap Mc becomes larger than the effective chemi- 
cal potential. We expect, that similar to the usual t — t' 
model the transition from a metal to insulator belongs to 
the universality class of Kosterlitz-Thouless transitionsi^ 

Let us now briefly comment on the behavior of the spin 
sector. Is is usefull first to consider the strong coupling 
limit U ^ A,t, t'. In this limit the low-energy physics 
of the t — t' ionic-Hubbard chain is described by the spin 
S" = 1/2 frustrated Heisenberg model 



eff 



— J y ] s„ • s„-i-i + J' y ^ s„ • s„ 



+2 : 



where the exchange couplings are given by 



J 



U 1-A2/C/2' 



(41) 



(42) 



For next-nearest neighbor couplings J' < 0.25J the spin 
excitation spectrum of the spin model (|1T|) is gapless and 
gapped for J' > 0.25 J.— Using (|42)) we easily conclude 
that at t' > 0.5t and f/ ^ t, A a spin gapped sponta- 
neously dimerized phase is realized in the ground state. 
Since for arbitrary finite alternating ionic potential the 



8 



ground state of the system is characterized by the pres- 
ence of a long-range ordered CDW pattern^S we conclude 
that the whole charge and spin gapped sector of the phase 
diagram at t' > 0.5< corresponds to a ferroelectric insu- 
lating phase. At t' < 0.5t and for strong repulsion the 
spin sector is gapless and therefore in this limit a Mott 
insulating phase is realized. Note, that since the ionic po- 
tential slightly enhances the exchange parameter J and 
does not influence (in first order with respect to t^ /U) 
the next-nearest-neighbor exchange J' for intermediate 
values of the on-site repulsion U > At, the Mott phase 
slightly penetrates into the t' > 0.5t sector of the phase 
diagram. 

In the weak-coupling limit at t' > 0.5t the FI phase 
undergoes a transition cither to the metallic phase or 
directly to the BI phase. Since the metallic phase with 
four Fermi points is characterized by a gapped spin sector 
for arbitrary weak on-site repulsion^i^^i^ we conclude 
that the spin gapped phase is a generic feature of the 
model for t' > 0.5i. For t' < O.bt, with increasing U the 
spin gap continuously decays and finally vanishes in the 
MI phase. 

To conclude our analysis we briefiy discuss the phases 
which are realized along the transition lines. The bor- 
der line between the BI and metallic phases corresponds 
to the Luttinger liquid state with gapless charge and 
spin excitation spectrum. The border line between the 
metallic and FI phases corresponds to the unconven- 
tional metallic phase with gapped spin and completely 
gapless charge excitation spectrum. The border line be- 
tween the BI and FI phases corresponds to the uncon- 
ventional metallic phase with gapless optical excitations 
and gapped spin and single-particle charge excitations. 
Finally, the border line between the FI and MI phases 



corresponds to the phase with gapped charge and gap- 
less spin excitation spectrum. 



IV. CONCLUSIONS 

We have studied the ground state phase diagram of the 
half-filled one-dimensional t — t' ionic-Hubbard model us- 
ing the continuum-limit bosonization approach. We have 
shown that the gross features of the ground state phase 
diagram and in particular the behavior of the charge sec- 
tor can be described by a quantum double- frequency sine- 
Gordon model with topological term. We have shown 
that with increasing on-site repulsion, for various val- 
ues of the parameter t', the model shows the following 
sequences of phase transitions: Band insulator - Fer- 
roelectric Insulator - Mott Insulator; Band Insulator - 
Nonmagnetic Metal - Ferroelectric Insulator and Non- 
magnetic Metal - Ferroelectric Insulator. 

We expect, that the transition sequence BI-metal-FI 
found in this paper is an intrinsic feature not only of 
the ID chain, but is a generic feature of the t — t' ionic- 
Hubbard model and will show up also in higher spatial 
dimensions. 

V. ACKNOWLEDGMENTS 

It is our pleasure to thank D. Bacriswyl, F. Gcbhard, 
D. Khomskii, R. Noack, A. Rosch and A.-M. Tremblay 
for many interesting discussions. GIJ and EMH acknowl- 
edge support from the research program of the SFB 608 
funded by the DFG. GIJ also acknowledges support from 
the STCU-grant N 3867. 



^ N. F. Mott, Metal-Insulator Transitions, 2"'' ed., Taylor 
and Francis, London (1990). 

^ F. Gebhard, The Mott Metal-Insulator Transition, 
Springer, Berlin (1997). 

^ M. Imada, A. Fujimori, and Y. Tokura, Rev. Mod. Phys. 
70, 1039 (1998). 

* For a recent review see F. H. L. Essler, H. Frahm, F. 
Gohmann, A. Kliimper, and V. E. Korepin, The One- 
Dimensional Hubbard Model, Cambridge Univ. Press, 
Cambridge (2005). 

^ A. Rosch, Preprint arXiv Condmat/0602656. 

^ F. Anfuso and A. Rosch, Preprint arXiv Cond- 
mat/0609051. 

A. Fuhrmann, D. Heilmann, and H. Monien, Phys. Rev. B 
73, 245118 (2006). 
^ M. Fabrizio, A. O. Gogolin, and A. A. Nersesyan, Phys. 
Rev. Lett 83, 2014 (1999); ibid, Nucl. Phys. B 580, 647 
(2000). 

^ A. P. Kampf, M. Sekania, G. I. Japaridze, and P. Brune, 

J. Phys. C 15, 5895 (2003). 
" S. R. Manmana, V. Meden, R. M. Noack, and K. 

Schonhammer, Phys. Rev. B 70, 155115 (2004). 



" H. Otsuka and M. Nakamura, Phys. Rev. B 71, 155105 
(2005). 

M. E. Torio, A. A. Aligia, G. I. Japaridze, and B. Normand, 
Phys. Rev. B 73, 115109 (2006). 

A. Garg, H. R. Krishnamurthy and M. Randeria, Phys. 
Rev. Lett. 97, 046403 (2006). 

P. Lombardo, R. Hayn and G. I. Japaridze, Phys. Rev. B 
74, 085116 (2006). 

S.S. Kancharla and E. Dagotto, Preprint ArXiv Cond- 
mat/0607568. 

M. E. Torio, A.A. Aligia, and H. A. Ceccatto, Phys. Rev. 
B 64, 121105(R) (2001). 
" A. A. Aligia, Phys. Rev. B 69, 041 101 (R) (2004). 

A. A. Aligia, K. Hallberg, B. Normand, and A. P. Kampf, 
Phys. Rev. Lett. 93, 076801 (2004). 

C. D. Batista and A. A. Aligia, Phys. Rev. Lett. 92, 246405 
(2004); tbid, Phys. Rev. B 71, 125110 (2005). 
^° N. Paris, K. Bouadim, F. Hebert, G.G. Batrouni, and R.T. 
Scalettar, Preprint arXiv Condmat/0607707 
J. Hubbard and J.B. Torrance, Phys. Rev. Lett. 47, 1750 
(1981); N. Nagaosa and J. Takimoto, J. Phys. Soc. Jpn. 
55, 2735 (1986): T. Egami, S. Ishihara, and M. Tachiki, 



9 



Science 261, 1307 (1993). 

V. J. Emery, in Highly Conducting One-Dimensional 
Solids, edited by J. T. De Vreese, R. P. Evrard and V. E. 
Van Doren (Plenum, New York, (1979); J. Solyom, Adv. 
Phys. 28, 201 (1979). 

A very exceptional case is the ID electron system with the 
short-range (SR) and long-range (LR) Coulomb interac- 
tion. As was shown in Ref. [23, in the case ol commensurate 
band-filling and SR interaction which support the forma- 
tion of a charge-density-wave (CDW) insulating phase, the 
increasing LR component of the Coulomb repulsion drives 
a transition from the CDW insulator to a metal. Further 
increase of the LR component leads to a transition from 
the metalic phase into the insulating Wigner crystal phase. 
D. Poilblanc, S. Yunoki, S. Maekawa and E. Dagotto, Phys. 
Rev. B 56, R1645 (1997). 

F. D. M. Haldane, Phys. Rev. B 25, 4925 (1982); ibid, 26, 
5257 (1982). 

K. Okamoto and K. Nomura, Phys. Lett. A 169, 422 
(1992); S. Eggert, Phys. Rev. B 54, 9612 (1996). 
S. R. White and I. Affleck, Phys. Rev. B 54, 9862 (1996). 
R. Resta and S. Sorella, Phys. Rev. Lett. 74, 4738 (1995); 
ibid, Phys. Rev. Lett. 82, 370 (1999). 

N. Gidopoulos, S. Sorella, and E. Tosatti, Eur. Phys. J. B 
14, 217 (2000). 

Y.Z. Zhang, C.Q. Wu, and H.Q. Lin, Phys. Rev. B 67, 
205109 (2003). 

M. Fabrizio, Phys. Rev. B 54, 10054 (1996). 

K. Kuroki, R. Arita, and H. Aoki, J. Phys. Soc. Japan 66, 

3371 (1997). 

S. Daul and R. M. Noack, Phys. Rev. B 58, 2635 (1998). 
R. Arita, K. Kuroki, H. Aoki, and M. Fabrizio, Phys. Rev. 
B 57, 10324 (1998). 

S. Daul and R. M. Noack, Phys. Rev. B 61, 1646 (2000). 
C. Aebischer, D. Baeriswyl, and R. M. Noack, Phys. Rev. 
Lett. 86, 468 (2001). 

K. Louis, J. V. Alvarez and C. Gros, Phys. Rev. B 64, 
113106 (2001); ibid, 65, 249903E (2002). 



^® K. Hamacher, C. Gros, and W. Wenzel, Phys. Rev. Lett. 

88, 217203 (2002). 
'■'^ M. E. Torio, A. A. Aligia, and H. A. Ceccatto, Phys. Rev. 

Lett. 67, 165102 (2003). 

C. Gros, K. Hamacher, and W. Wenzel, Europhys. Lett. 
69, 616 (2005). 

''^ M. Capello, F. Becca, M. Fabrizio, S. Sorella, and E. 
Tosatti, Phys. Rev. Lett. 94, 026406 (2005). 
G.I. Japaridze, R.M. Noack and D. Baeriswyl, unpub- 
lished, condmat/0607054. 

L. Balents and M. P. A. Fisher, Phys. Rev. B 53, 12133 

(1996) . 

A. O. Gogolin, A. A. Nersesyan, and A. M. Tsvelik, 
Bosonization and Strongly Correlated Systems, Cambridge 
Univ. Press, Cambridge (1998). 

T. Giamarchi, Quantum Physics in One Dimension, 
Clarendon Press, Oxford (2004). 
■'^ G. I. Japaridze and A. A. Nersesyan, JETF Pis'ma 27, 356 
(1978); [JETP Lett. 27, 334 (1978)]; ibidem, J. Low Temp. 
Phys. 37, 95 (1979); V. L. Pokrovsky and A. L. Talapov, 
Phys. Rev. Lett. 42, 65 (1979); H. J. Schuiz, Phys. Rev. B 
22, 5274 (1980). 

G. I. Japaridze, A. A. Nersesyan, and P. B. Wiegmann, 
Nucl. Phys. B 230, 511 (1984). 

Al. B. Zamolodchikov, Int. Jour. Mod. Phys. A 10, 1125- 
1150 (1995). 

""^ S. Lukyanov, A. Zamolodchikov, Nucl.Phys. B 493, 571 

(1997) . 

^° R.F. Dashen, B. Hasslacher and A. Neveu, Phys. Rev. D 
10 3424 (1975); A. Takhtadjan and L.D. Faddeev, Sov. 
Theor. Math. Phys., 25, 147 (1975). 

G. Delfino and G. Mussardo, Nucl. Phys. B 516, 675 

(1998) . 

^2 F. D. M. Haldane, Phys. Rev. B 25, 4925 (1982); K. 
Okamoto and K. Nomura, Phys. Lett. A 169, 433 (1993); 
G. Castilla, S. Chakravarty, and V. J. Emery, Phys. Rev. 
Lett. 75, 1823 (1995). 



